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A charged parallel plate capacitor will create particle-antiparticle pairs by the Schwinger process and 
discharge over time. We consider the full quantum discharge process in 1+1 dimensions including 
backreaction, when the electric field interacts with massless charged fermions. We recover oscillatory 
features in the electric field observed in a semiclassical analysis and find that the amplitude of the 
oscillations falls off as t -1 / 2 and that stronger coupling implies slower decay. Remarkably, Ohm's 
law applies to the vacuum and we evaluate the quantum electrical conductivity of the vacuum to 
be 2e/y/Tr, where e is the fermionic charge. Similarities and differences with black hole evaporation 
are mentioned. 



I. INTRODUCTION 

Following Schwingcr's 1951 paper [l[, it is well-known 
that quantum effects cause electric fields to produce 
electron-positron pairs. The effect can be interpreted as 
the tunneling of virtual electron-positron pairs into real 
particles. One imagines a virtual e + e~ pair produced in 
the vacuum which is then torn apart by the background 
electric field, say within a capacitor, with the positron 
accelerated in the direction of the electric field, and the 
electron in the opposite direction. A similar effect has 
been proposed in de Sitter space where the rapid ex- 
pansion of spactimc "pulls" particles out of the vacuum. 
Hawking radiation from black holes has also been inter- 
preted Schwinger process, though the universally 
attractive nature of gravity suggests that there are dif- 
ferences. For example, in the electromagnetic case, it 
is clear that once the charges are created, the positive 
charge accelerates away from the positively charged ca- 
pacitor plate due to electromagnetic repulsion. In the 
black hole case, however, there is radiation even though 
the black hole attracts all outgoing particles. 

The energy for pair creation in an electric field must 
eventually come from the energy in the electric field itself. 
Hence the electric field has to decay due to the Schwinger 
process, just as Hawking radiation is assumed to cause 
black holes to evaporate. The problem of electric field de- 
cay clearly involves calculation of the backreaction of the 
Schwinger process and this is a hard problem. There have 
been several attempts to analyze the decay of the elec- 
tric field by semiclassical methods, replacing quantum 
operators by their expectation values [2| ■ The results are 
interesting. For example, a uniform electric field will not 
discharge monotonically but will undergo oscillations. If 
the conclusion can be directly transported to the black 
hole case, it would imply oscillations of the black hole 
mass and not monotonic evaporation. A key difference 
though is that electric charges can be positive or nega- 
tive, whereas the particles in Hawking radiation can only 
have positive mass. 

In this paper we re-visit the problem of capacitor dis- 
charge due to the Schwinger process, without restrict- 



ing ourselves to the semiclassical approximation. We can 
solve the full quantum problem but the price we pay is 
that we are then only able to treat massless fermions 
and the exponential suppression of the classic Schwinger 
process is absent. 

We treat the case of massless QED in 1+1 dimensions 
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The fermions interact with the gauge field by the stan- 
dard minimal coupling, and an electric field leads to 
fermion pair production. The advantage of 1+1D QED 
is that it can be bosonized to yield , 
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where, e M „ is the Levi-Civita tensor in 1+1D with eoi = 1 
and 
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where, without any loss of generality, we can assume 
g > 0. The correspondence between the fermionic and 
bosonic models is given by_ identifying the currents at the 
quantum operator level 
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The bosonized model is particularly easy to solve be- 
cause it is quadratic in fields and hence there are no in- 
teractions. It is sufficient to solve it classically. This 
simplification only occurs if the fermions are massless in 
the original model. If we had included a mass for the 
fermions, we would have obtained a sine- Gordon poten- 
tial for the scalar field, which is an interacting scalar field 
theory, and a quantum treatment of the bosonized model 
would become necessary. 

To analyze the discharge of a capacitor, we would 
like to set up an initial electric field that is localized 
within a finite region of space and then examine its evolu- 
tion due to the spontaneous production of fermion pairs. 
This leads to a physical difficulty because the capacitor 
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plates are necessarily external to the system. This phys- 
ical difficulty arises even in classical electromagnetism, 
where Maxwell's equations are solved but the boundary 
conditions are provided externally. The difficulty can 
be avoided in gravitational systems, for example, dur- 
ing gravitational collapse, because the system naturally 
tends to evolve toward a black hole. On the other hand, 
if the gravitational problem is set up to include an eternal 
black hole as an "external" agent, a similar issue arises. 
There is no analog of gravitational collapse in the elec- 
tromagnetic situation since similar electric charges repel. 
If we were to set up a configuration of <f> field correspond- 
ing to separated localized positive and negative charges, 
the charge distributions would simply spread out due to 
mutual repulsion and then annihilate. Hence an external 
capacitor plate is necessary to set up the problem. 

Here we will set up the capacitor problem in two ways. 
In the first scheme, we fix the boundary condition that 
there is a charge +Q at x = — L/2 and another charge 
— Q at x = +L/2. These external charges are taken to 
be fixed and are non-dynamical. The Schwinger process 
then creates pairs and dissipates the electric field within 
the capacitor but cannot "evaporate" the charge on the 
plates. In this case, we find that the capacitor charges 
are screened due to the Schwinger process and the electric 
field decays exponentially with distance from a capacitor 
plate. Our second strategy to set up the capacitor plates 
is to introduce an external potential (two "bags") into 
which we can insert the charges +Q and — Q as config- 
urations of the 4> field itself. These charges are now dy- 
namical because the field cj> is dynamical. The Schwinger 
process causes evaporation of the charge from the bags. 
In both physical realizations, the approach to the asymp- 
totic state is not monotonic but oscillatory. In the first 
case, the final static state still contains the ±Q charges 
but the charges are screened by opposite charges. In the 
second case, the final state is not static. Instead it con- 
tains bound states of positive and negative charge densi- 
ties that oscillate without dissipation within the bags. 

In the following two sections we describe the discharge 
of a capacitor in the two physical setups, first with "ex- 
ternal charges" and second with "external potential" . We 
can solve the first setup analytically, allowing us to ob- 
tain explicit expressions for the late time behavior of the 
current, electric field, and energy decay law. The so- 
lution of the second setup has only been obtained nu- 
merically. The results of both methods are summarized 
in Sec. IIVI and show that the capacitor discharge is os- 
cillatory, the root-mean-square current is proportional 
to the root-mean-square electric field (Ohm's law), and 
that the electrical conductivity of the vacuum is equal to 
g = 2e/y/ir. 



II. SETUP I: EXTERNAL CHARGES 
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FIG. 1: Schematic view of setup in section UT] Two infinitely 
heavy charges, +Q and —Q are placed at x = —L/2 and 
x = +L/2 respectively. The initial electric field is given by 
the thick dark line. The expected final electric field is given 
by the dashed grey line. 



tively (see Fig.Q]). The electric field due to these charges 
satisfies Maxwell's equations and is a non-zero constant 
in the region between the plates. The value of the elec- 
tric field in the region —L/2 < x < +L/2 is E = Q. For 
| a; | > L/2, the electric field vanishes. This electric field 
configuration, together with <fi — and <j) = dt<t> = 0, 
corresponding to no particles, are the initial conditions 
whose evolution we will consider. 

The equations of motion follow from ([2]). Including the 
external charges on the capacitor plates, we get 
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with 



f; = ge^d^ (7) 
j» xt = Qu v {5{x-L/2)-5{x + L/2)) (8) 

and u v = (1,0). 

Maxwell's equations © can be integrated immediately 
to get 



F 



where 



F = Q(@(x + L/2) - Q(x - L/2)) 



(9) 



(10) 



Note that we have set the constant of integration to zero 
so that the electric field at spatial infinity vanishes. 

The charge within some interval (a, b) is given by 
Gauss' law 



In this section we treat the capacitor as made up of 
two external charges ±Q placed at x = ^L/2 respec- 



F i (t, x = b) - F Q1 (t, x = a) 

g[<f>(t,x = b)-<t>(t,x = a)] (11) 
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Inserting ([9]) into the scalar field equation, ([5]), gives 

(d 2 + m 2 )^ = -mF 

= -mQ(e(x + L/2) - Q(x - L/2)) (12) 

where the effective mass of the scalar field is given by the 
coupling constant, 



m = g 



(13) 



Hence our problem reduces to solving the Klein- Gordon 
equation for a scalar field of mass to, sourced by the 
"electric field", E\ — i*oi = Q, within the capacitor. 
The initial condition at t = is given by the requirement 
that no fermions be present, or in terms of the bosonic 
variables, 



4>{t = 0,x) =0 = <j>(t = 0,x) 



(14) 



Before we evolve the equations, however, it is inter- 
esting to find static solutions into which the system can 
evolve asymptotically. 



A. Static Solution 

In the asymptotic future, t — > oo, we expect the <j) 
solution to be simply the static solution to (IT21 . Since 
the Klein-Gordon equation (fl"2"|) is linear, we may first 
solve it with F = (Q/2)(&(x) - Q(-x)) = (Q/2)sgn(x) 
i.e. due to a single point charge at the origin. The static 
solution to the present problem would then follow using 
appropriate linear superposition. For now, 



(d 2 +m 2 )<j> 
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Using the integral representation of the step- function and 
that of the retarded Green function, G r (x — y), 



Gr(x-y) 



Q(x) 
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a ik-(x— y) 



(2tt) 2 -{k Q - t0+) 2 + (k,) 2 + m 
provides us with the integral representation 

<f>(x) = -to| / d 2 y G T (x - y)(Q(y) - e(-y)) 



(17) 



-mQ 



dk sin(fcx) 
27r k(k 2 + to 2 ) 



(18) 



The integral in (|18p may be evaluated by perform- 
ing a partial fraction decomposition of the denomina- 
tor k(k 2 + to 2 ) and converting the resulting three sub- 
integrals into appropriate contour integrals, which may 
then be computed straightforwardly. The answer is 



cj) (x) 



2m 



sgn(x) ( 1 
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The static solution to (|T2|) is therefore 

4> s {x) = ~Mx~L/2) + M* + L/2) 

_Q / er™\*\ sinh(TOL/2), \x\ > L/2 
to \ 1 - e- mL / 2 cosh(mx), \x\ < L/2 K > 

B. Dynamical Solution 

The solution to (IT21 we are seeking must satisfy the 
initial conditions in (1141) . To obtain this dynamical 
solution, we add a homogeneous solution, </>h, obeying 
(d 2 + TO 2 )0h = 0, to the static solution <fi s such that the 
initial conditions are satisfied. Again, it helps to first 
solve the problem with a single charge. Then we have 
to solve (|15|) with the initial conditions corresponding to 
(fl"4|) . From the conditions in (fT4|) . we observe that the 
integral representation of the solution is 



(j){t,x) 



<f) s (x) + (f)h(t,x) 

dk sin(fca;) 
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1 -c 
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The solution to (|12|) with the same initial conditions is 
thus 



<t>(t, x) 



= -2m 



<t>(t, x - L/2) + (f>(t, x + L/2) 
dk cos(fcx) sin(fcL/2) 
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From this, we can extract the flux of energy passing 
through a given spatial point x and integrate over all 
time to get the total energy radiated. It is 



/•OO 

F(x) = dt T c 
= -(2mQ) 2 



dk_ 

27 



dt do4>di(f> 
dp sin(fcL/2) cos(kx) 

2^ k(k 2 + TO 2 ) 



shi(pL/2) sin (pa;) 



We first use 

r +co dk exp(ikx) _ 
J-oo 2^ k 2 -a 2 ~ 

followed by 

/dk exp(ikx) 1 
2^fc 2 (fc 2 + a 2 ) = ~2^ 

for a > 0, to obtain for \x\ > L/2 



1 sin(a|o;|) 

2 a ' 
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and, for |x| < L/2, 
O 2 

T(x) = ^- (2mx - e~ mL sinh(2mx)) . (27) 
4m 

By considering the limit \x\ — > oo, we obtain the total 
energy that is radiated 



T la .d = 2T{oo) 



Q 2 L 



1 - — (1 - e~ mL ) 
mL v ' 



(28) 



We can check that the expressions for F{x) are consistent 
with energy-momentum conservation, d x T 01 = — 6\T 00 , 
if the final field configuration is the static <p s . A direct 
calculation verifies 



d x T{x) = T 00 {t = 0, x) - T Q °(t = oo, x) (29) 



with 



T m = \idtcj,) 2 + kd^) 2 + \ (mcf> + F) 2 (30) 



where the "mass term" in T 00 arises from the electric 
field energy density (1/2)Fq 1 , since Foi = n%4> + F. 

One may also use formula 3.876.1 in Gradshteyn and 
Ryzhik 



dk- 



sm 



\/k 2 + a 2 



■ cos(fefc) 



fJo^vV— r/j, 0<b<p 
0, b>p>0 
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and apply it to the integral representation (|22[) , to obtain 
the rate of charge creation, dtj° oc dtd x (j>: 



d t d x <t>{t,x) = — 



J [mJt 2 -x 2 _ 0(t- 



- J [mJt 2 -x 2 + )Q(t-x+) 



(32) 



where x± = x ± L/2. The first term can be attributed 
to the —Q charge at x = +L/2 whereas the second to 
the +Q charge at x = —L/2. Charge creation at any 
location x goes to zero at late times because the Bessel 
functions tend to Jo(m\t\) and hence dtj° —> 0. This is 
consistent with the expectation that the asymptotic field 
configuration is <p s (x) in ([2H)l . 

The electric current at any spacetime location is 

■x 1 r. 2^ f dk cos(fcx) sin(fci/2) . 
j = gc/) = -2m z Q / — — '- — sm(u k t) 
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_ (33) 
where we have introduced uik = Vk 2 + m 2 and also used 
g = m. We use a trick to evaluate this integral. Let us 
first differentiate with respect to I = L/2. Then after ap- 
plying appropriate trignometric identities and Eq. ()31l) . 



this gives 
dif -- 



2m 2 Q 



„ cos(kx) cos(fcZ) . , 
dk i — — sm(uj k t) 
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+ J (rn^Jt 2 -x 2 _^j Q(t - x+) 
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where we have defined x± = x ± I. Noting that the 
current vanishes when the plate separation vanishes (I = 
0), we get 



f{t,x) = f L/ ' 



J mJt 2 -x 2 _ )Q(t-x-) 



+ Jo (rnyjt 2 -x%j Q(t - x+) 



(35) 



At late times, t 3> \x±\, we can Taylor expand the 
Bessel functions at mt and then perform the integration 
over I to get 
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(36) 

where we have used the asymptotic form of the Bessel 
function Q in the second line. The first term is a good 
approximation for t 3> mL 2 , \x±\, m -1 . 

The expression in Eq. (|36[) shows that the current 
within the capacitor (say at x = 0) oscillates at the "mi- 
croscopic" frequency given by m. If we average out these 
fast oscillations, the cosine gets replaced by 1/V2 and we 
find that the root-mean-squared current decays as t^ 1 ! 2 : 
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The electric field within the capacitor decays to a 
static value that can be obtained from the static solution 
Eq. (f2"U| inserted into ©. The time-dependent electric 
field within the capacitor at late times can be obtained 
from the expression for the current in Eq. p6|) together 
with the asymptotic static solution 



E{t, X) = -Estatic 
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To check this expression simply differentiate with respect 
to time and keep the leading order term in l/{mt). This 
agrees with E = g<fi = j x . Note that the static part of 
the electric field plays no role. In fact, well away from 
the capacitor plates, the static electric field dies off expo- 
nentially fast and can be ignored. So we will define the 
decaying part of the electric field as E^ = E — E static and 
refer to this as the electric field. 



The expression for the electric field shows that it is 90° 
out of phase with the current but the amplitude has the 
same decay as the current. The rms value of the 

electric field is 



E, 



d.rms 



This leads to Ohm's law 
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where cte is the electrical conductivity of the vacuum 
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(T E = 9 



(41) 



This result is independent of Q and L. 

Let us now consider the energy in the capacitor. At 
late times, the fields approach the static solution whose 
energy can be computed using ([20)1 
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T oo dx =9-h- e - 2Lr 
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This shows that the final energy is smaller if the coupling 
g = m is stronger. Or stronger coupling implies more 
complete radiation of the capacitor energy. 

We can identify a typical time scale for energy loss from 
the capacitor by considering the ratio of the decaying 
part of the energy in the electric field within the capacitor 
at time t to the initial energy (Q 2 L/2). The ratio is 



£(t) m 2 Q 2 L 3 /(8irmt) _ 



8(0) 



where 



Q 2 L/2 
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(43) 



(44) 



Hence larger couplings imply longer decay times i.e. 
slower decay. The capacitor is more effectively discharged 
when the coupling constant is large but it takes a longer 
time for the discharge to happen. In the zero coupling 
limit, the rate of pair production is rapid, but the original 
electric field E remains relatively undissipated. 



III. SETUP II: EXTERNAL PLATES 

In the second setup we do not wish to introduce ex- 
ternal charges. Instead the capacitor is charged with the 
same fermionic field, i/j, or its bosonized version, <f>. How- 
ever we still need to have some capacitor "plates" that 
we can charge. These plates have to be external to the 
system. To implement this scheme, we add a double well 
potential to the action in © 



S V = -^J d 2 x[V(x + L) + V(x - L)]cj) 2 



(45) 




FIG. 2: Schematic view of setup II. The dashed line is the 
double well potential V(x + L) + V(x — L). The thick black 
line is the initial <j> configuration, tanh(x- + L) — tanh(x — L), 
and the thick grey line is a bound state solution in the double 
well. Note that the figure is meant to be schematic and the 
horizontal axis is not the zero of the potential. 



where the form of V(x) is chosen so that we can find 
a non-dissipative solution for the scalar field in a single 
well. A sketch of the setup is shown in Fig. [2] 
It is convenient to choose 



V(x) = -2M 2 sech 2 (Mx) 



(46) 



where M is some mass scale that sets the depth and width 
of the well. With this choice the solution to the single 
well problem 



{d 2 + m 2 + V(x))<f> = 
is given by the Xcj> 4 kink 

<t>{t,x) = cos(mi) tanri(Mx), 



(47) 



(48) 



where we have chosen the initial condition <j){t = 0, x) = 
0. The solution describes a positive charge in the well at 
t = 0, which then oscillates - due to pair creation in the 
electric field - but does not dissipate. The charge in the 
well at any given time can be found from Eq. © to be 
2m cos(mi). 

Now consider a capacitor with two plates, one at x = L 
and the other at x = —L. Since we do not have any ex- 
ternal charges, we take the constant background electric 
field, F in Eq. (|9]), to be zero. This means we are now 
solving 

(d 2 +m 2 + V(x + L) + V(x - £))</> = (49) 

with the following choice of initial conditions 

<t>(t = 0,x) = tanh(M(a; + L)) - tanh(M (x - L)) = <f> (x) 
4>(t = 0,x) = (50) 

In what follows, we shall set M = 1 and so all quantities 
will be in units of M. 
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A. Asymptotic State 

Before solving for the time evolution, we consider the 
asymptotic state, which will be a stationary solution of 
Eq. (|4T))) . That is, we think of the double well equation 
(|49|) as a Schrodinger equation 



(51) 



with Hamiltonian 



H = -d 2 + m 2 + V(x + L) + V(x - L) (52) 

The corresponding Hamiltonian for a single potential 
well centered at x = is 



H x = -fig + m 2 -+ 
This has exactly one bound state 



V(x) 



tp oc sech(x) 



(53) 



(54) 



with eigenvalue ui 2 — m 2 — 1 Q. (Recall that we are 
working in units with M = 1.) So for the double well 
potential, at least when the two wells are well separated, 
there must be exactly two bound states, which can be 
approximated as 



ip ~ sech(x — L) ± sech(x + L) 



(55) 



The energies of these two bound states are nearly iden- 
tical, uj± = m 2 — 1 ± e~ r , T{m,L) ^> 1, split by ex- 
ponentially small corrections due to tunneling between 
the double wells. These bound states are the stationary 
states that the system can evolve into. 

We also observe that there is an apparent instability 
in the current model when m < 1, since then luq < 
and the bound state solution can grow exponentially. To 
understand this instability, we examine the double well 
action in (|45|) . If the potential V(x) is deep enough, there 
will be a region where m 2 + V(x) is sufficiently negative, 
that it becomes favorable for <fi to grow without bound 
in this region. In terms of the fermionic model, the well 
is so deep that it is favorable to pull fermion pairs out of 
the vacuum and put them at the bottom of the potential. 

The evolution of the initial data, 4>q, in this setup can 
be evolved formally by writing 



(j)(t,x) 



COS (Uj n t) l/jn{x)(ljjn\^o) 



(56) 



where the summation is over both bound and continuum 
states of H . We can check that the initial conditions are 
satisfied by setting t = in the factor cos(w n i) and in its 
time derivative. We expect that, as time progresses, the 
continuum states will disperse to infinity, leaving behind 
only the initial overlap with the bound state. While for- 
mally correct, the expansion in Eq. (|56[) is only useful if 
we know the full eigenspectrum of the double well poten- 
tial. In the absence of the eigenspectrum, it is easier to 
numerically evolve the equation of motion. 
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FIG. 3: Envelopes of the plots of the current at the center 
of the capacitor, j x (t,x = 0), versus time. From black to 
light grey, the curves represent, respectively, the evolution 
for m — 1.5, 2 and 2.5. The rapid oscillations between the 
envelopes are not shown. 
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FIG. 4: Envelopes of the plots of the electric field at the 
center of the capacitor, E(t, x = 0) versus time. From black 
to light grey, the curves represent, respectively, the evolution 
for m — 1.5, 2 and 2.5. The rapid oscillations between the 
envelopes are not shown. 



Time Evolution 



We have evolved Eq. (|49|) using the explicit Crank- 
Nicholson algorithm with two iterations with first-order 
absorbing boundary conditions. The runs were done on 
very large lattices so that boundary effects are minimal. 

In Fig. [3] we plot the current envelopes at x = versus 
time for several different parameters, disregarding the 
rapid oscillations between the envelopes. Similarly, in 
Fig.2]we show the behavior of the electric field at x = 0. 
On the log-log plot in Fig. [SI it is clear that the envelopes 
decay as a power law. A fit gives 



'A 



(57) 



where TV sa 14 is a factor which could depend on the 
dimensionless product ML, where M _1 is the width of 
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FIG. 5: Envelopes for the currents (solid curves) and the elec- 
tric fields (dashed curves) at x = on a log- log plot. Com- 
parison with the light straight line with slope —1/2 clearly 
shows the l/\ft fall off. From black to light grey, the curves 
represent, respectively, the evolution for m = 1.5, 2 and 2.5. 
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FIG. 6: Ratio of the envelopes of current (divided by m) to 
electric field at x = for m = 1.5, 2 and 2.5. The constant 
flat ratio implies Ohm's law and the value of 1 implies that 
the conductivity is g. 



the wells (see Eq. (HH)). The electrical conductivity is 
therefore again given by Eq. (|41 [) . as for setup 1. 



IV. CONCLUSIONS AND DISCUSSION 

Wc have studied the quantum discharge of a capacitor 
in masslcss QED in 1+1 dimensions by bosonizing the 
model. The bosonized model is non-interacting and can 
be solved classically. The solution includes all backreac- 
tion effects. We now summarize the key results. 

The final state depends on the setup used to describe 
the capacitor plates. We have chosen two different ways 
to describe the capacitor plates. In both cases, the fi- 
nal state is non-trivial. In setup I, the plates keep their 
original charge but the charges are screened due to the 
Schwingcr process. In setup II, there are no external 
charges, but there are external potentials that play the 



role of capacitor plates. Then the final state consists of 
fermion-antifermion pairs that are bound to the plates. 

The energy in the final state depends on the coupling 
constant, g, and equivalently the mass of the scalar field, 
to. The final state energy decreases with increasing g, 
while the time for the capacitor to discharge increases 
with increasing g, as seen in Eq. (|44"]) . So stronger cou- 
pling leads to more complete discharge but the discharge 
process itself is slower. We suggest that the longer dis- 
charge time for larger coupling constant is due to the 
tighter binding of fermion-antifermion pairs that need to 
be split apart by the electric field. 

The discharge process is highly oscillatory, as also seen 
in the semiclassical analysis [2[ and the amplitude of os- 
cillations falls off rather slowly, as i -1 / 2 . This suggests 
that the massless QED system is underdamped. 

Our results show that the root-mean-square current in 
the capacitor is directly proportional to the root-mean- 
square electric field, indicating that Ohm's law holds on 
a macroscopic scale. Thus it makes sense to define the 
electrical conductivity for the massless QED vacuum to 
be oe = jxmsl ^rms and our results indicate the simple 
relation cte = g which can also be written in terms of the 
fermionic charge as cte = 2e/y / 7r. 

A correspondence is often made between Schwinger 
pair creation and Hawking radiation, though we have 
indicated differences between the two processes that 
prompt us to use caution in drawing a correspondence. 
If the oscillatory features of the discharge process carry 
over to black hole evaporation, we may expect black hole 
mass oscillations during evaporation. Though, in con- 
trast to the capacitor, the black hole system is unstable 
in that smaller mass black holes are hotter and evapo- 
rate faster, while weaker electric fields in the capacitor 
do not discharge faster. So it would appear that a fluc- 
tuation that excessively decreases the mass of the black 
hole, would make it evaporate yet faster in what may be 
a runaway process. 

The issues of black hole formation and the final state 
of black hole evaporation cannot be resolved by this cor- 
respondence since the capacitor plates have to be intro- 
duced externally, whereas there are no such externally 
set conditions in the case of gravitational collapse. Yet 
it would be extremely interesting if the electromagnetic 
Ohm's law has a gravitational analog that relates energy 
flow (current) from a black hole, or during gravitational 
collapse, to the "gravitational electric" field (see Sec. 4.4 
of Q). Perhaps the instability of the black hole can be 
summarized in a negative "gravitational conductivity" . 

A potential application of our findings is to supercon- 
ducting cosmic strings, where massless QED in 1+1 di- 
mensions is expected to apply for fermion zero modes on 
the string [l(| • Our analysis shows that if superconduct- 
ing strings really behave as 1+1 dimensional systems, 
they will carry oscillatory currents because of the back- 
reaction of the induced currents on the external electric 
fields. (Oscillatory currents were also discussed in 
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though these occurred due to the periodic dynamics of 
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the strings.) 

While our analysis has enabled us to fully treat back- 
reaction of the Schwinger process, our results cannot be 
transported to 3+1 QED for two reasons. First, the elec- 
tron has a non-zero mass. For electric fields smaller than 
the electron mass squared, the Schwinger process is ex- 
ponentially suppressed and the vacuum is essentially an 
insulator. In situtations where the electric field is larger 
than the electron mass squared, the exponential suppres- 
sion is absent and the dynamics may be closer to what we 
have found. The second reason is that the larger number 
of dimensions can change the picture dramatically. In 
1+1 dimensions, the inter-charge potential is linear and 
electric charge is confined. In 3+1 dimensions, electric 
charges interact by the Coulomb potential and are not 
confined. This suggests that our system may be closer to 
the case of chromo-electric fields in 3+1 dimensions with 
massless quarks. While this has some features that re- 
semble the model we have considered, there are essential 



differences due to the non-Abclian nature of the interac- 
tions. 
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